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INTRODUCTION 



A complete understanding of the many puzzles related to black holes, such as entropy 



will require a quantum theory 



[lj], Hawking radiation [2], and the final state question 
describing gravitational collapse, where both matter and geometry are fully dynamical. The 
ideal would be the ability to follow the evolution of an initial matter-geometry quantum state 
to "black hole formation" and beyond. This goal has been implicit in many works on black 
hole physics beginning with Unruh thirty years ago followed a few years later by work 
of Q. More — it n as ta ^ as one of the motivations 

for studying quantum gravity [5|. To date, however, no complete quantum framework is 
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available for studying gravitational collapse of a scalar field. 

This paper continues a recent line of work aimed at developing the tools necessary to 
realize this ideal. So far, we have developed the quantum kinematics for spherically sym- 
metric gravitational systems, in a setting which leads to a result concerning the resolution 
of the classical singularity at the quantum level j^. Furthermore, we have proposed how to 
define a black hole at the quantum level, without fixing any classical "horizon boundaries" 
jjj. Finally, we formulated the classical dynamics of black holes in a framework that results 
in a true Hamiltonian, rather than a system with a Hamiltonian constraint P| . 

The current paper is concerned with the development of the quantum dynamics of the 
gravity- scalar field system in spherical symmetry. It consists of two main parts. The first 
one deals with the conceptual issues connected with the notion of true dynamics. This is 
necessitated by the fact that our treatment of the dynamics is based neither on Dirac's 
constraint quantization, nor on a complete gauge-fixing, but lies in between those two ex- 
tremes by employing a partial gauge-fixing only of time. As this type of approach has - to 
the best of our knowledge - not been discussed before in any detail, we include a careful 
discussion. The second part is of a more technical nature and describes the construction of 
the Hamiltonian operator starting from the reduced HamUtonian derived in Q. 

The work is an attempt to complete a proposal for the quantum theory of the gravity- 
scalar field system in spherical symmetry. It provides a calculational framework for studying 
gravitational collapse in quantum gravity, which we hope will lead to a better understanding 
of Hawking radiation, and of the final state of gravitational collapse. 

The next section discusses some conceptual issues concerning quantum evolution in a 
setting where only a partial (time) gauge fixing is utilized. This is followed in Section III by 
a review of the kinematical framework for quantization of the 1+1 dimensional field theory 
that describes the gravity-scalar field system. Section IV contains a construction of the 
Hamiltonian operator. Since the classical Hamiltonian contains a square root, we introduce 
here a new method of defining the corresponding operator. In Section V we give the action 
of the Hamiltonian on basis states, describe how to implement unitary evolution, and discuss 
the issue of dynamical singularity resolution. We conclude in Section V with a summary 
and outlook for numerical implementation. 
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II. CONCEPTUAL SETTING 

There are at least two methods to introduce a notion of time into a theory with time 
reparameterization invariance. The most obvious way is by means of gauge-fixing, where a 
suitable function on the classical phase space is taken to be the time function. Another is 
the closely related method of using partial and complete observables where a degree 



of freedom is chosen as the reference clock, and the evolution of the remaining degrees of 
freedom is measured with respect to that clock. 

We have chosen the former since it is closest to the setting in which the well-known 
semiclassical results about black holes have been derived. In earlier work jj| on the classical 
theory, we derived a reduced Hamiltonian for the gravity-scalar field theory in spherical 
symmetry by performing only a time gauge fixing, and leaving the remaining spatial coor- 
dinate freedom untouched. The Hamiltonian constraint was solved as a strong condition, 
and a reduced spatial diffeomorphism generator remained as the only (first class) constraint. 
Together with the surface term, this constraint forms the total Hamiltonian. If the surface 
term is written as a bulk integral and combined with the other bulk terms, one can identify 
the true local Hamiltonian density and a diffeomorphism generator for the remaining degrees 
of freedom. 

Fixing only a time gauge raises the question of the extent to which unambiguous evolu- 
tion can be achieved in the quantum dynamics, since evolution still contains a gauge part 
manifested in the freedom in the shift function. This may be seen schematically through 
the requirement that a gauge condition such as f(p, q) = be preserved in time. It leads to 
the equation 

= / = {/,/ (NH + N a H a )}, (1) 

which fixes the lapse function N in terms of N a . The latter remains arbitrary, and must be 
specified to compute evolution classically or quantum mechanically. This situation is similar 
to Yang-Mills theory where the total Hamiltonian density before gauge fixing and solving 
the Gauss law is 

H YM = ~{E 2 + B 2 )+A i G i , (2) 

where A* is the Lagrange multiplier and G % is the Gauss law expression. 

In contrast the Hamiltonian for our problem (in a fixed time gauge) contains an additional 
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twist (see below for the details). Schematically, it is of the form 

H grav = f ((N r )'; q, vr) + N r C r (q, vr), (3) 

where C r is the radial diffeomorphism generator, q, 7r denote the collection of canonical phase 
space variables, and N r is a lagrange multiplier (the remaining radial component of the shift 
function in spherical symmetry). The similarity between the two cases (YM and gravity) is 
the clear separation of the gauge generators from the "true" Hamiltonians. The difference 
is that the gravity case contains Lagrange multiplier dependence also in the first term, but 
only through its radial derivative. 

This poses the conceptual issue of obtaining unambiguous evolution, since the gravity 
reduced Hamiltonian generates a family of time evolutions for observers specified by N r . 
This freedom is limited to the spatial reference system only. Once it is fixed so is the clock. 
This situation lies in the middle of the two common scenarios in dealing with Hamiltonian 
gravity, where either no gauges are fixed, or all are fully fixed. 

We will address this issue in two steps. Firstly, we obtain an Hamiltonian operator that 
depends on N r , determine its action on basis states, and thereby obtain a prescription for 
evolution with any N r (with the prescribed asymptotic conditions). This gives a family of 
evolved states parameterized by N r . A unique evolution is then obtained by specifying a 
specific functional form for this function. This may be viewed as the quantum analog of 
evolving with a fixed shift function. 

Secondly, we will not impose the diffeomorphism generator as a constraint on states 
a'la Dirac. Rather we will give prescriptions for obtaining information pertaining to the 
collapse problem in a manner which is manifestly invariant under this symmetry. This 
means looking at the dynamics of certain phase space observables, such as null expansion, 
as functions of other phase space variables. Suitable semiclassical states that are peaked 
on classical solutions of the radial diffeormorphism constraint will be used as the physical 
states. 

Taken together, the implementation of these ideas, provided in the following sections, 
appears to provide a tractable approach to the gravitational collapse problem in quantum 
gravity. 
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III. REVIEW OF THE KINEMATIC AL SETTING 



Before addressing the issue of quantum dynamics, we briefly review the quantum kine- 
matical setting, as introduced in p|. The starting point is a classical field theory in 1+1 
dimensions, characterized by the canonical pairs (R,Pr), describing the geometry degrees 
of freedom, and (0, P^), describing a scalar field. The basic variables that are turned into 
quantum operators are the smeared fields 



CO 



Rf = I Rf dr (4) 

and 



o 



CO 



dr, (5) 



where / and g are suitable test functions. These configuration variables, together with the 
translation generators 

U x (P R )(r) = e^M (6) 

and 

U x {P^){r) = e iXp ^\ (7) 

form a closed Poisson algebra, and are taken as the basic variables to be converted into 
operators. All other operators are be constructed in terms of these. 

The Hilbert space for the quantum theory is spanned by basis states of the form 

= \a x . . . a Nl ]b x . . . b N2 ), (8) 

where cik,bi are real numbers, and N% and N 2 are positive integers. The factors of L in 
the exponents reflect the length dimensions of the respective field variables. The intuitive 
picture is that each basis state 'tests' the quantum scalar field at iV points in space and the 
basis of the 'excitation space' at each point consists of (generalized) plane waves. The inner 
product is 

(ai . . . ajvi; bi . . . bN 2 \a'i ■ ■ ■ a Nl \ b x . . . b' N2 ) 

= $ai,a' 1 --- S b N2 ,b> ' , (9) 
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if the states contain the same number of sampled points, and is zero otherwise. The action 
of the basic operators is 

Rf | oi . . . cin-l ; b\ . . . b N2 ) = 
^^akfixk)]^ ...a Nl ;bi. ..b Na ) 



(fig \ai . . . a Nl ; b\ . . . bjv 2 ) — 
I?^big{yi)\ax . . . a Nl ;bx . . . b N% ) 



(10) 



and 



e i ^ p ^)\a 1 ...a Nl ;b 1 ...b N2 ) = 
\ai . . . , aj — Xj, . . . ajvj ; b\ . . . b^ 2 ) 
e VkP*(tt,)| 0l . . . a Nl ; bx . . . b N2 ) = 

\a\ . . . ajv^ bi . . . ,bk — fik, ■ ■ ■ &iv 2 ) 



(11) 



where aj (resp. bk) is if the point Xj (resp. yk) is not part of the original basis state. In 
this case the action creates a new 'excitation' at the point Xj (resp. y^) with 'mode' — Aj 
(resp. — fik). These definitions give the commutator 



R f , e iAP «M 



■A/(x)L 2 e lAP «W. 



(12) 



Comparing with the classical Poisson bracket relation, and using the Poisson bracket com- 
mutator correspondence, it turns out that L = a/2/p, where lp is the Planck length. A 
similar commutator relation holds for the matter scalar field. 

In this formalism we can also define other operators of interest. For example for the 
Hamiltonian, we need operator analogs of the radial derivatives R' and 0', and the square 
of the momentum P|. 

Definitions for the operators corresponding to R' and other derivatives of fields are ob- 
tained by implementing the idea of finite differencing using the operator Rf We use 
narrow Gaussians with variance proportional to the Planck scale, peaked at coordinate 
points rk + elp, where < e 1 is a parameter designed to sample neighbouring points J^: 



fe{r, r k ) 



1 



exp 



[r - r k 



eh 



21% 



(13) 
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Denoting Rf c by R € for this class of test functions we define 

R'{r k ):=^- (A £ -Ro). (14) 

A further motivation of this form is that in the gauge R = r the corresponding classical 
expression gives unity in the limit e — * 0. This definition captures the simplest finite 
difference approximation to the derivative at the operator level. Second derivatives may be 
similarly defined by converting a finite differencing scheme into the corresponding operator. 

The quantization defined above does not give definitions for operators corresponding to 
momenta, since only translation generators are directly realized as operators at the first 
step. This of course is similar to any quantum theory on a lattice. Operators for momenta 
can however be realized using the translation generators, for example by expressions such 
as 



and 



" 'U X -U{\ (15) 



^ 2 :=| (2-*Wl) (16) 

These A dependent expressions will be utilized below for the Hamiltonian operator. 

Finally, let us note that the background independent (ie. metric free) quantization out- 
lined here is not the same as the one for the scalar field reported in |l2 . ll^ . where the 
basic variables used are the integral of the scalar field momentum over space J s P^, and the 
exponential of the scalar configuration e lX ^. The present approach has the advantage that it 
is fairly straightforward to write a local Hamiltonian density operator using the translation 
and <f)' operators. A further technical point of difference, apart from the choice of basic 
variables, arises from the fact that R and <fi are scalars, and Pr and P^ are scalar densities 
of weight 1. This means that the functions / in (jlj) and the A in (JHJ) are scalar densities 
of weights 1 and —1 respectively. These combine to give eig envalues of Rf and 0/ that are 
scalars. In contrast, for the quantization studied in |l2|, [l3(, it is the space-integrated mo- 
mentum density that is diagonal; this is the direct analog of the surface integrated densitized 
dreibein in loop quantum gravity. 

This ends our review of the kinematical setup. Our goal in the next section is to define 
the Hamiltonian operator for the collapse problem on this Hilbert space. 
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IV. HAMILTONIAN 

The classical Hamiltonian in a time gauge fixing corresponding to flat spatial slicing 
given was derived in The classical phase space before gauge fixing has an extra pair of 
canonical variables (A, Pa), with the spatial metric 

ds 2 = A 2 dr 2 + R 2 dn 2 . (17) 

The condition A = 1 is second class with the Hamiltonian constraint, which is solved clas- 
sically for the conjugate momentum Pa- This leads to the reduced Hamiltonian 

POO 

Hi = / [(N r yP A + N r (P R R' + P^')}dr 
Jo 

= J (N r )' (rP r + ^{RP R f-X) dr 

POO 

+ / N r {P R R' + P^') dr, (18) 
Jo 

where N r is the remaining (and still arbitrary) non-zero component of the shift function A^ a 
after imposing spherical symmetry, 

X = 16R 2 (2RR" - 1 + R' 2 ) + 1QR 2 H^ (19) 

and 

* = # + (20) 

The Hamiltonian (fTH|) is obtained by writing the surface term in the reduced action as a 
bulk integral and combining terms. It gives the time-gauge fixed evolution equations for the 
fields R(r,t) and <p(r,t), and their canonical conjugates. Our main goal in this paper is to 
construct the corresponding operator. 

Let us write the first part of the Hamiltonian density (excluding the radial diffeomorphism 
term) as 

H^ = h 1 + ^h 2 -h 2 -h 2 + n = hi + A (21) 

where we have defined 

h = P R R } K = -1QR 2 (2RR" - 1 + R' 2 ) (22) 
h\ = 8P|, h\ = SPY 2 - (23) 
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There are two ways to construct an operator corresponding to a square root. One ap- 
proach is to find the eigenvalues of its argument and work with the corresponding basis 
of eigenvectors. The operator can then be defined as the square root of the eigenvalues if 
the spectrum is positive semi-definite. A potential alternative is to see if the square root 
operator can be defined using Dirac's idea. 

In this procedure an operator for the classical function H = ^/p 2 + m 2 (i = 1 • • • 3) 
is constructed by writing H = aipi + (3m using anticommuting matrices /3, such that 
H 2 = (pi + m 2 )I, where H 2 defined by matrix multiplication. This works because the 
momentum components pi commute. This is not true for the elements hi in the argument 
of our Hamiltonian, so Dirac's trick cannot be used, at least in its basic form. 

There is the additional problem that the term 7Z, which is the Ricci scalar, is not a 
squared quantity, and so can be negative. (In the gauge R = r this term vanishes since 
we are using the flat slice asymptotic conditions [9].) H is however diagonal in the basis 
since the operator analogs of the fields R, R' and R" are all diagonal 8]. Therefore at least 
for the subset of states where its eigenvalues X-n are positive we can define an operator h/± 
whose eigenvalues are \fX%. We will return to this issue below after describing an approach 
to address the square root problem. For now, and the following discussion, we define 

h\ := U. (24) 

The idea is to obtain a definition of the square root part A of the Hamiltonian operator 
by working in a larger Hilbert space 

n = n Kin ® v, (25) 

for some V to be specified, and writing 

A = h k e\ (26) 

where the h\ act in Tixin, and the e k act in V. The key requirement is that the action of A 
in Tt must be such that 

A 2 = h k k V kl I, (27) 

where r] kl = diag(H h) and I is the identity operator in the space V. All other operators 

O acting in Hxin are extended to TC by the identity action in V. Since the hk do not 
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commute, this means that the operators e k must satisfy 

e k e l = T] kl I. (28) 

That no such operators exist for any space V may be seen by the following argument. Let 
\ijj k ) be a complete set of normalisable states, and let us assume that operators satisfying 
(|28j) exist. Then we are led to contradictions such as 

i = (m = w {e 1 ) 2 W 

= (^(e 4 ) 2 ^) = We 1 g*X)l^>^*l g * gl lV'> = 0- (29) 

k 

There is a way out of this situation if we demand the weaker condition that 

(^|i 2 |^> = ^\h k k r) kl \V), (30) 

for every |$) G H that is of a form specified below. The basic idea is to introduce a 
four- dimensional Euclidean vector space V and consider its decomposition into orthogonal 
subspaces 

V = ® k V k (31) 

k — 1 • ■ • 4. Let us denote by P k the projection operators onto these subspaces. By definition 
the Pfc satisfy P k Pi = Pk^ki- Now define operators e k by 

e k = P k (32) 



for k = 1,4, and 



for k = 2, 3. These give 



e k = iP k (33) 



A 2 = h k hie k ei = h k hiP k r] H (34) 
Consider now the action of A 2 on states of the form 

|*) = |V)|P> (35) 

where e Uk™, 

|P>=5>*>' ( 36 ) 
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and the \pk) denote the basis of V corresponding to the decomposition into the The 
result is 

A2 w = (mm) pm) 

k,i 

= J2(h h hk\iI>))\Pk). (37) 
k 

From this it is evident that ()3()j) holds. Finally, we still need a refinement to get a self-adjoint 
A. To do this we define the operators h% and h^ such that they have the property 

h\ = h 2 h\ = h 3 . (38) 

These ensure that the square root operator A in the extended Hilbert space is self-adjoint. 

The last issue to address is the definition of the operator corresponding to 1Z (the Ricci 
scalar) in Eqn. (|22p. This contains derivatives of R which are defined by the finite difference 
operators introduced in jsj. Since these are diagonal on basis states, so is 7Z . Now, in 
the gauge R = r, which fixes the shift vector N r to be proportional to l/v 7 ^ everywhere, 
the Ricci scalar vanishes. Our Hamiltonian is not in this coordinate gauge (which is still 
free), but as alluded to in Sec. II, the idea is to evolve quantum states with this choice of 
lapse function, with initial states for which the values of the radial field R are distributed 
linearly with the graph points r^, ie. ~ rv This guarantees that the eigenvalue of 7Z 
vanishes, at least initially. However it will not remain so because the Hamiltonian contains 
terms with Pr, which is represented by the operator 

P R (r k ) := ^ (e ap ^ - e~ ap ^) . (39) 

The translation operators on the r.h.s. act to shift excitations at the point which result 
in the eigenvalue of TZ being moved from zero at the selected point. However this move 
is very small since < A << 1. Because of this we define the action of by 

km = v\>^m, (40) 

where is a basis state. This definition makes the assumption that, starting from a basis 
state with X-ji = 0, the deviation from zero after the action of Pr is positive. 

This completes our prescription for the Hamiltonian operator for the scalar field collapse 
problem. With each of its constituents well defined, it is straightforward to compute its 
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action on a basis state. The fact that this can be done in contrast to the case for full gauge 
fixing represents some progress which we anticipate will lead to concrete calculations for 
quantum collapse. A strategy to do this is presented in the next section. 

We close with a discussion of the quantization choices, or ambiguities, inherent in the 
steps outlined here. Firstly, let us note that the general procedure used in obtaining the 
square root operator is new. The only other possible approach is to seek the spectrum of the 
argument of the square root, which given its form appears a formidable task. If this could 
be done, it would represent another, possibly physically distinct, choice of Hamiltonian 
for this problem. This circumstance would be analogous to the Dirac and Klein-Gordon 
Hamiltonians. Secondly, h\ contains products of non-commuting operators so a choice must 
be made for its operator version. The natural one is to take the symmetric product (RPr + 
PrR)/2. Finally, the only other ambiguity in the definition of the Hamiltonian is in the 
choice of "lattice" parameter A. It is natural to also use the same parameter in the Gaussian 
smearing functions used in operators such as Rf. In implementing evolution numerically, the 
goal of course is to ensure that physical results do not depend on this parameter. Like any 
numerical computation, one would like to see that the evolution of a fixed initial state leads 
to a convergent answer for the final state after a fixed number of time steps. This means 
that when the computation is repeated for successively smaller values of A, the answers for 
evolved physical variables have asymptotic "continuum" values. 



V. QUANTUM EVOLUTION 

In classical numerical simulations of scalar field collapse the general approach is to start 
with initial data represen ting a shell of scalar field, and to evolve it using some choice of 
lapse and shift functions JJJ. At each step of the simulation an "apparent horizon" check 
is made as a criterion for black hole formation. This is a null geodesic trapping condition at 
each step or leaf of the evolution. In spherical symmetry the goal is to find the outermost 
radial location on each leaf where the condition is satisfied. The evolution of this location 
is taken to represent the dynamical boundary or horizon of the black hole. 

There are at least two concrete versions of what is a dynamical horizon. The first was 
formulated by Hayward Q|. In addition to the usual criteria for null geodesic expansions, 
this work consists of additional conditions designed to distinguish future, past, inner and 
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outer local horizons. The second lifts some of these conditions, and points out that the 
resulting definition of dynamical horizon allows a nice formulation of local flux laws. The 
common and minimal feature of both definitions are the conditions 

6+ = 0, < (41) 

where 6± are the in(out)going null geodesic expansions. Equivalent information is captured 
in the observable # + #_, which goes from negative to positive as a dynamical black hole 
boundary is crossed. 

With the Hamiltonian operator defined in the last section, and the operator analogs of 
the null expansion operators given in [8|, we are in a position to give a procedure for a 
quantum collapse calculation. This involves specifying (i) initial states, (ii) an evolution 
procedure, and (iii) a quantum test for black hole formation. 

The first question is what are suitable initial states. The basis states represent values 
of the scalar field 0(r, t) and the radial field R(r, t) at a set of discrete coordinate points 
7*1 • ■ -Tk, which is a sample of the half line. As such these states may be compared with 
the discrete data for a classical numerical simulation. This suggests the use of "profile" 
states where we take for example the scalar field excitations to have a gaussian profile, and 
the radia 
coherent 

These would be infinite linear combinations of basis states, and so computation with them 
would be more involved, especially in the present field theory setting. 

The second step is the implementation of evolution. Rather than constructing a finite 
evolution operator by exponentiation of the Hamiltonian, which is very cumbersome, it is 
more suitable to implement repeated infinitesimal evolution using a suitable scheme. The 
simplest possibility 

W)t+At = (I-iAtHm t (42) 

is not unitary. However there are unitary schemes available for this purpose. One example 
is based on the well known Crank-Nicholson method, where the Schrodinger equation in 
discrete time (labelled by n) is written as 

^ (|VWi - |V>„) = \ (H\iP) n+1 + %> n ) . (43) 



excitations to be linearly distributed. Another possibility is the use of suitable 
is| or other form of semi-classical states that are peaked at classical configurations. 
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This leads to the manifestly unitary (but implicit) evolution scheme given by 

1 + |V>n+i =(l-^AtH) |V)n (44) 

Since H depends also on the free (classical) function N r , so does the evolved state \ip) n+ \. 
This function must be fixed (with the fall off condition N r ~ 1/ yjr) to get a unique evolution 
jjj. The simplest choice is to take this form for all points r/~ in the chosen initial state. 

The third step is to implement a test for black hole formation at each time step of 
the evolution. As mentioned above, the minimum requirement for this is that we must 
have operators corresponding to the null expansions. A prescription for constructing these 
were given in ^]. In addition we require this test to be invariant under the remaining 
radial diffeomorphism constraint. This is achieved, for example, by looking at the quantities 
(9 + (rk,t)9-(rk,t)) as functions of < i2(r&,t) >, where the expectation values are in the 
state arrived at by stepwise evolution from some initial state. The resulting curve is radial 
diffeomorphism invariant, and its intersection with the (R) axis gives the location and size 
of the evolving horizon. It is the dynamics of this curve which is of interest for black hole 
formation and subsequent evolution at the quantum level. 

There are many other useful diffeomorphism invariant quantities of interest that can be 
computed at each time step. Two examples are the scalar field configuration (<p(r k ,t)), and 
a curvature measure such as (fi(r k ,t)) (the trace of the ADM momentum) j]| , both viewed 
as functions of (R(ri-,t)). 

The second and third steps are to be repeated for multiple time steps to extract time 
dependent profiles of the functions of interest, such as the ones just mentioned. It is in 
this manner that physical information about collapse may be obtained in the kinematical 
Hilbert space after fixing the time gauge classically. Although evolution is with respect to a 
fixed N r so that the evolved states depend on it, the physical information contained in the 
suggested functions is independent of the coordinate points r k : radial diffeormorphisms act 
to shift the points r k to r' k without changing the values of the field variables. 

VI. DYNAMICAL AVOIDANCE OF THE SINGULARITY 

In earlier work we gave a construction of an operator corresponding to the classical 
variable 1/R that is bounded on the Hilbert space we are using for the quantum theory . 
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This result has the direct consequence that curvature singularities in spherical symmetry 
are avoided. This is because phase space variables that classically diverge at the singularity 
do so as positive powers of 1/R. The corresponding quantum operators are constructed 
as products of the 1/R operators, and so are also bounded. This does not involve the 
Hamiltonian in any way so the result may be viewed as kinematical singularity avoidance. 
In this section we show, using the construction of the Hamiltonian operator given above, 
that inclusion of dynamics does not alter this result. 

The fundamental question here is whether quantum evolution remains well-defined 
through the region of highest curvature, or whether it stalls or breaks down there. Clas- 
sically dynamics is encoded either in (i) the Hamiltonian constraint, or (ii) a Hamiltonian 
derived from a partial (time) gauge-fixing, or (iii) a Hamiltonian derived from a complete 
gauge-fixing as in ^j]. Although our work is concerned with the second case, we will look 
at the issue of singularity avoidance from all three points of view. Let us consider first the 

Hamiltonian constraint. Its classical expression is 

1 



H 



R 2 A 



1(PaA) 2 -\(PaA)(PrR) 



+— \2RR"A - 2RR'A' - A 3 + AR' 2 ] 

A 2 L J 



+ 



P 2 



R 2 

+ 7TT 



J'2 



2AR 2 2A Y ^ 
Due to the 1/R factors, it is divergent at the points r where R — 0. Upon quantization 

these factors turn into bounded operators as shown in []|. From this result, and the form 
of the basic quantum operators it is clear that the quantum Hamiltonian constraint has 
finite action on all states, including those for which the eigenvalue of Rf is zero. The action 
of the constraint operator on a state with maximum eigenvalue of the 1/R operator gives 
a linear combination of basis states that contain shifts in excitation values of the fields 
generated by the action of the corresponding momentum operators (J39|) . It is not difficult 
to see that due to this, each of the states in the linear combination correspond to a lower 
eigenvalue of the 1/R operator. Thus in this sense a time step "evolution" by the action of 
the Hamiltonian constraint on a state of maximum value of curvature gives a new state in 
which its expectation value is lower than the maximum. 

Consider next the fully reduced Hamiltonian that was found in Q] 

r j. 



dr 



4r : 



2 j/2 



exp 



S<j,(r')dr' 



(46) 
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where 

This is clearly divergent at r = 0. As r is a parameter rather than a configuration field 
variable, quantization of the totally reduced theory cannot resolve the singularity at the 
quantum level. 

Finally, let us consider the Hamiltonian whose quantization is the subject of this paper. 
Firstly, an inspection of the classical reduced Hamiltonian (|18|) shows that, unlike the Hamil- 
tonian constraint, it has no manifestly divergent 1/R factors. This surprising feature is just 
a consequence of the time gauge fixing, and concomitant strong solution of the Hamiltonian 
constraint. (If one continues the reduction process further by the gauge choice R = r, the 
divergence reappears via the shift function, which is proportional to 1/y/r This becomes 
the source of the divergence at r = in the fully reduced Hamiltonian.) Secondly, the 
Hamiltonian operator has well defined action on basis states, and evolution does not stall on 
states with maximum eigenvalue of the 1/R operator. Rather, it gives a linear combination 
of states each of which has a lower eigenvalue of this operator (for the same reason as for the 
action of the Hamiltonian constraint discussed above). In closing this section, we compare 
the above scenario for dynamical singularity resolution with the cosmological case studied 



m 



3,0 



and the Schwarzschild black hole studied in |2l| . In both these cases the systems 
are finite dimensional, unlike the model in this paper. In the first case, the action of the 
Hamiltonian constraint on a basis state gives a finite difference equation with coefficients 
such that action on the state of zero volume gives a bounce. In the second case the interior 
(Kantowksi-Sachs) and exterior of the Schwarzschild spacetimes are quantized with appro- 
priate matching at the event horizon, which is taken as the fundamental classical dividing 
line. Our perspective is that these cases are quite different from the matter coupled case 
we treat here. In particular, with non-vanishing matter fields, the extended Schwarzschild 
spacetime does not arise, so quantizing it to discuss singularity resolution is not relevant for 
the quantum treatment of the collapse problem. 

VII. DISCUSSION 

We have constructed the quantum Hamiltonian for the spherically symmetric system of 



gravity coupled to a scalar field in a fixed time gauge. Together with earlier papers 
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this work completes the construction of a quantum theory for studying the gravitational 
collapse of a scalar field in spherical symmetry. We are now in a position to address the 
questions surrounding black holes that have been generally acknowledged to find a resolution 
only within a full quantum treatment. 

Any application of our formalism to a given physical situation will require a choice of 
initial state. As all the black holes that have been detected so far have macroscopic size, 
the first task is to determine how to represent what we know as a classical black hole in the 
quantum theory. This obviously calls for a construction of the semiclassical sector, including 
he search for suitable semiclassical states. This will be the subject of a future publication 
2^ |. Once the issue of what initial state to take has been settled, one can then investigate 
the quantum time evolution of the system. 

Among the problems of physical interest is the evolution of a quantum state that satisfies 
the quantum horizon conditions. One would like to know how this evolution depends on 
the matter part of the state, and how the horizon grows or shrinks. An indispensable tool 

n 

for such questions will be the null expansion operators which were constructed in [qj, which 
serve as horizon finders. It should be pointed out that, even if one does not subscribe to 
the exact definition of the horizon of a black hole proposed in |8|, these operators would 
still play an important role in any other approach to a quantum description of black hole 
horizons. 

Perhaps the most interesting questions to address in this framework are whether and how 
Tawking radiation arises, and what is the end point of collapse. A preliminary investigation 
~v\ . indicates that the end result of gravitational collapse is a Planck size remnant. This 
appears to be intimately connected with the existence of an upper bound on curvature, 
which in turn has an appealing intuitive analogy with Fermi pressure; whereas the latter can 
be ultimately overcome by gravitational forces, the former cannot since it is a fundamental 
quantum gravity effect. It represents the ultimate limit to which matter can be compactified. 
Finally, as time evolution in our setting is unitary, one can surmise already at this point 
that the solution to the so-called information loss paradox is that there has never been a 
paradox in the first place. 

Apart from applications, another potential line of investigation is to improve on the 
framework developed here. One of the questions here concerns time gauge fixing. There are 
of course many other choices, so it would be useful to see if there are others that might lead 
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to simpler reduced Hamiltonians. One possibility is to look at only the time gauge fixing 
used in jj], without fixing the radial gauge R = r. Another approach to the same issues 
using the connection-triad variables is being developed by Bojowald and collaborators [23]. 
This work offers a parallel approach in the loop quantum gravity programme. 
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